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Abstract: An approach is presented which allows a self-consistent description of 
the fragmentation of single-particle strength for nucleons in finite nuclei employing 
the Greens function formalism. The self-energy to be considered in the Dyson 
equation for the single-particle Greens function contains all terms of first (Hartree- 
Fock) and second order in the residual interaction. It is demonstrated that the 
fragmentation of the single-particle strength originating from the terms of second 
order can efficiently be described in terms of the so-called BAGEL approximation. 
Employing this approximation the self-energy can be evaluated in a self- consistent 
way, i.e. the correlations contained in the Greens function are taken into account 
for the evaluation of the self-energy. As an example this scheme is applied to 
the nucleus 16 O, using a realistic nucleon nucleon interaction. The effects of the 
correlations on the occupation probabilities and the binding energy are evaluated. 
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Most of the microscopic nuclear structure studies are based on the independent particle 
model (IPM) describing the nucleus as a system of nucleons moving without any cor- 
relations in a mean field (Hartree-Fock potential) generated by the interaction with all 
other nucleons. A lot of effort has been made, both form the experimental as well as the 
theoretical side, to explore the borderlines of this simple shell-model and search for clear 
fingerprints of correlations beyond the IPM. Nucleon knock-out experiments with electrons 
or hadrons have been performed [1-4] to determine absolute spectroscopic factors and study 
deviations of occupation probabilities from the predictions of the IPM. 

The depletion of occupation numbers for hole-states, states which are completely 
occupied in the IPM, partly originates in the strong short-range components of a realistic 
nucleon-nucleon (NN) interaction. The effects of these short-range correlations have mainly 
been studied for nuclear matter [5-8] assuming that the main results can be transferred 
to finite nuclei. Here we want to focus the attention on long-range correlations. This 
implies that one has to consider a finite nucleus since such correlations are sensitive to the 
low-energy excitation spectrum, reflecting the shell-structure of finite nuclei. The effects 
of short-range correlations are taken into account by evaluating an effective interaction 
appropriate for the model-space under consideration. As described in ref.[9] this is done 
by solving the Bethe-Goldstone equation for a realistic OBE potential[10] assuming a Pauli 
operator appropriate for this model-space. 

As a starting point of the discussion let us consider the definition of the irreducible 
self-energy for the nucleons taking into account the terms of first and second order in the 
effective interaction V e //. The first order contribution can be written 

E S = E/ ^ < «7|Ve//|/V > G^M , (1) 

77' 

while the terms of second order are defined by 




X < ail\Veff\l'5' X 1/6\VeffW > Gjyiu-U! + U 2 )G S S' (uJt)G ^ (u 2 ) .(2) 

The summations on single-particle quantum numbers 7,5... are restricted to orbits which 
define the model-space. In our test calculation on 16 (9 discussed below, single-particle 
states up to the lpOf shell are taken into account. In equations (1) and (2) G a p(uj) refers 
to the Greens function which is obtained from the Dyson equation 

G a p(uj) = S a pg a (uj) +^2g ce (uj)T,^(u)G 1 f 3 (uj) (3) 

7 
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and g a is the Greens function determined for the self-energy contribution of first order 
in the interaction E^ 1 ). It is obvious that equations (1) to (3) have to be solved in a 
self-consistent way. As a first step of an iteration scheme to reach this self-consistency we 
consider the Hartree-Fock (HF) approximation for g a 
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It is diagonal in the basis of HF single-particle states and defined in terms of the single- 
particle energies e^ F for states a above and below the Fermi energy e F . In the first iteration 
step also the Greens functions G in eqs.(l) and (2) are replaced by the corresponding 
HF Greens functions. On this level the self-energy terms E^ and E^ 2 ) are represented 
by the Goldstone diagrams of figure la) and lb) plus lc), respectively. Inserting this 
approximation for the self-energy into the Dyson equation one can obtain the solution by 
solving the following^ set of equations (see ref.[9,ll,12]) 
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where all elements of the matrix that are not indicated, are zero, 
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3 "Pi 1 -P2 

is the energy of a 2 particle 1 hole (2plh) configuration j and 

dj =< pip 2 \V\ah > , 



(6) 



(7) 



the coupling of the single-particle state a to the 2plh configuration j. The normalization 
of the eigenvectors in eq.(5) to unity gives the correct normalization of the spectroscopic 
factors (Xq ) 2 . In a similar way the quantities E^ and bk refer to 2 hole 1 particle (2hlp) 
configurations. The Greens function G a is then given in the Lehmann representation as 
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(8) 



t In order to keep the notation simple we ignore here and in the following presentation 
the non-diagonal terms in the self-energy. 
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Due to the coupling of the single-particle state a to the 2plh and 2hlp configurations the 
single-particle strength is distributed to various discrete states at energies u n . In general 
one will observe for single-particle orbits a, which in the IPM are below the Fermi level, 
an occupation probability 

n a = ^Te(e F -c^) W) 2 , (9) 

n 

which is different from one. Similarly, one obtains a non-zero occupation probability for 
single-particle states above the Fermi level. 

In a second iteration step one may try to employ the Greens function of eq.(8) in the 
definition of the self-energy. This is no problem for the £W. Rather than the conventional 
BHF single-particle energies we will now obtain 

e% HF =< a\t\a > + ^ < a(3\V eff \a(3 > np . (10) 

P 

Translated into the language of diagrams: The modification of the Greens function due to 
the second order terms in the self-energy yields a depletion of the occupation of the hole 
states, as represented in lowest order by diagram Id), and an occupation larger than zero 
for "particle" states as expressed by le). This correction is similar to the Renormalized 
Brueckner Hartree Fock approach ( therefore we use the label RHF in eq.(10)) accounting 
for depletion of holes and occupation of particle states[8,13]. 

It requires much more effort to use the improved Greens function also in the definition 
of the second order self-energy (see as examples the diagrams of figures If) and lg)). Again 
the solution of the Dyson equation (3) can be transformed into a problem of a matrix diag- 
onalisation. However, replacing the single-pole approximation of the HF Greens function 
of eq.(4) by the Greens function of eq.(8) leads to an explosion of basis configurations in the 
matrix equation which will replace eq.(5) in this second iteration step. Each combination 
of energies u n for the various orbits a yields a new configuration to be considered in this 
matrix. This increase of the number of configurations reflects the fact that in this second 
iteration step beside 2plh configurations also 3p2h, 4p3h etc configurations are taken into 
account as indicated by the diagrams If) and lg). 

Van Neck et.al. have proposed a scheme to represent the distribution of the single- 
particle strength in terms of strength selected in energy bins[14]. Here we would like to 
follow a different route, based on the assumption that for practical purposes the Greens 
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function of eq.(8) is well approximated by a small number of poles, which are representative 
for the spectral distribution at low energy. To introduce the method we recall that eq.(5) 
may be rewritten in a basis in which the matrix is tri-diagonal 
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This tri-diagonal form of the matrix is obtained by applying the so-called BAGEL (BAsis 
GEnerated by Lanczos method[9,15]) procedure to single-particle or single-hole states \i. 
In eqs (5) and (11) K and L denote the total numbers of 2plh and 2hlp configurations 
which have the same symmetry quantum numbers (parity, isospin and angular momentum) 
as the single-particle orbit ji. In this BAGEL basis we may truncate the matrix displayed 
in eq. (11). We call it the BAGEL(p,q) approximation if, in addition to the 1-body state \i 1 
we restrict the basis to include only the first p combinations of 2plh configurations (with 
diagonal elements Q\ . . . Q p ) together with the first q combinations of 2hlp configurations 
(Oi . . . Og). If \x refers to a hole state in the IPM we will consider p combinations of 
2hlp configurations plus q combinations of 2plh configurations in addition to the single- 
particle state \i. Diagonalization in this truncated basis yields p+q+1 eigenvalues and 
we may define the Greens function in the BAGEL(p,q) approximation by restricting the 
summation in eq. (8) to the eigenvalues generated in the truncated basis. 

It is obvious that BAGEL(0,0) yields the HF approximation for the Greens function 
and BAGEL(iiT,L) refers to the untruncated set of equations (5). The BAGEL(0,1) ap- 
proach is the simplest approximation leading for each fj, to an eigenvalue below and above 
ep which consequently yields non-trivial results for the occupation probabilities defined in 
eq.(9). It is easy to verify that the BAGEL(p,p) approximation reproduces the moments 
of the spectral distribution 

S m = J2<( X of ^ (12) 



evaluated for the Greens function of the untruncated equations from order m = up to 
the order m = 2p + 1. Therefore this BAGEL scheme provides a set of approximations 
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for the Greens function in terms of p+q+1 poles, which can systematically be improved 
ranging from the HF approximation (1 pole) to the result of the untruncated eq.(5). As 
an example figure 2 displays the distribution of single-particle strength (Xq ) 2 obtained for 
the Osl/2 state in 16 (9 in various BAGEL(p,q) approximations. 

The key point of the procedure presented here is that a BAGEL(p,q) approximation 
is inserted in the calculation of the self-energy (2). This self-energy then contains a much 
larger number of poles than when the single-pole approximation(4) is inserted in (2). The 
Dyson equation may still bee written in a matrix form (5). The dimension of the matrix is 
much larger. Therefore an exact solution is out of the question. The BAGEL truncation 
scheme can still be applied. Truncating again with the same values for p and q, an iterative 
scheme is obtained, which yields a self-consistent solution of eqs.(l-3). As an example let 
us consider the single-particle state a = Osl/2 for w O in the model-space discussed above. 
In the first iteration step, i.e. approximating G by the HF or BAGEL(0,0) approach the 
dimension is K+L+l = 82. If now for the second iteration step we employ the BAGEL(p,q) 
approximation for (p,q) = (0,l),(l,l),(2,2) or (6,6) the corresponding dimension is raised to 
596, 2178, 9638 or 160558, respectively. But even for a matrix of dimension 160558 it 
is not very difficult to determine the BAGEL(6,6) approximation to be used in the next 
iteration step and continue until convergence is achieved. 

In order to demonstrate the sensitivity of the results to the choice of p and q in 
this scheme, table 1 exhibits some typical results obtained in a self-consistent calcula- 
tion employing various BAGEL(p,q) approximations. For the hole states we present the 
occupation probabilities n a as defined in eq.(9) and the mean value for the hole-energy: 



while for selected particle states we show the corresponding mean value for the particle 
energies. Furthermore also the calculated binding energy per nucleon 



The correlation effects included in this scheme produce remarkable differences com- 
pared to the HF (or IPM) approach. The occupation probabilities for the hole states are 
reduced by about 3 percent and the single-particle energies for these states are more at- 
tractive by around 1 to 2 MeV. The total binding energy is increased by around 1 MeV per 
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is given. 
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nucleon. A fast convergence is obtained comparing the different BAGEL approximations. 
Already the simplest approach beyond the IPM, the BAGEL(0,1) scheme, yields a result 
very close to the most sophisticated one (p=q=6) , which we have been studying. 

In order to demonstrate the importance of self-consistency, the sixth column of table 
1 shows results obtained in the BAGEL(6,6) approach after the first iteration step, i.e 
considering the HF approximation for the Greens function in calculating the self-energy. 
This implies that only terms as displayed in figures la) - lc) are taken into account. One 
finds that this first iteration tends to overestimate the importance of correlations. 

The last column of table 1 shows the single-particle energies obtained in the RHF 
approximation according to eq.(10). At this level diagrams of figure la), Id) and le) are 
taken into account. One finds that the inclusion of the occupation probabilities according 
to Id) and le) increases the single-particle energies of all states as compared to HF. The 
additional inclusion of diagrams displayed in figures lb), lc), If) and lg) decreases the 
energies for the hole states but increases the energies of particle states. Therefore one 
enhances the gap between mean values for particle and hole energies, as defined in eq.(13), 
at the Fermi level. This should not be confused with the fact [16] that correlations decrease 
the gap in the excitation spectrum, as those energies are related to the solutions of the 
Dyson equation that are closest to the Fermi energy. 

The results of these studies demonstrate that the BAGEL scheme introduced here 
yields a very powerful tool for the representation of the single-particle Greens function, 
which can be used to solve the problem of calculating self-energy and Greens function in 
a self-consistent way. This scheme should produce reliable results also for heavier nuclei 
and model spaces larger than the one considered here. Larger model spaces and inclusion 
of residual interaction between the 2plh and 2hlp configurations tend to enhance the 
correlation effects [9]. A scheme very similar to this approach for finite systems could also 
be developed for studies of infinite systems like nuclear matter. 
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Table 1: 

Results for the single-particle energies (see eq.(13)) and the occupation probabilities for 
selected states in 16 as obtained in self-consistent BAGEL(p,q) calculations. Also listed 
are the energies per nucleon calculated according to eq.(14). The last 2 columns contain 
the results obtained in BAGEL(6,6) after the first iteration and the RHF single-particle 
energies according to eq.(10). All energies are listed in MeV. 
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Figure Captions 

Figure 1: Contributions to the self-energy (a, d, e) and E^ 2 ) (b, c, f, g). While the diagrams 
a) - c) show the contributions obtained in the first iteration step (replacing the Greens 
functions in eq.(l) and (2) by the HF approximation), the diagrams d) - g) show 
examples of contributions generated through the self-consistent scheme. 

Figure 2: Distribution of the single-particle strength {{Xq ) for the various poles u n in the 
Greens function) as obtained in different BAGEL(p,q) approximations for the 0si/ 2 
state in 16 0. Note the logarithmic scale and the fact that the strength obtained in 
BAGEL (1,1) at positive energies is hidden by the corresponding result obtained in 
BAGEL(0,1) as both are essentially identical. 



